Nonsensitive nonlinear homotopy approach 
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Generally, natural scientific problems are so complicated that one has to establish some effective 
perturbation or nonperturbation theories with respect to some associated ideal models. In this 
Letter, a new theory that combines perturbation and nonperturbation is constructed. An artificial 
nonlinear homotopy parameter plays the role of a perturbation parameter, while other artificial 
nonlinear parameters, of which the original problems are independent, introduced in the nonlinear 
homotopy models are nonperturbatively determined by means of a principle minimal sensitivity. 
The method is demonstrated through several quantum anharmonic oscillators and a non-hermitian 
parity-time symmetric Hamiltonian system. In fact, the framework of the theory is rather general 
that can be applied to a broad range of natural phenomena. Possible applications to condensed 
matter physics, matter wave systems, and nonlinear optics are briefly discussed. 
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In general, natural scientific problems cannot be ex- 
actly described by some analytical expressions. What 
scientists usually do is just first to establish some ideal 
models, and then to consider additional effects by use of 
some perturbation and nonperturbation theories. 

Perturbation theories (PTs) [l], Q already have quite a 
long and interesting history in physics. They have been 
used to explore physical systems that can not be solved 
exactly but with a small parameter. Though great chal- 
lenges are coming from the blooming numerical calcula- 
tion methods, PTs are still showing their great power 
in solving various problems. For instance, perturba- 
tion theories have been used to study the spectrum of 
non-hermitian parity-time {VT) symmetric Hamiltonian 
0, the structure of trapped Bose-Einstein condensates 
(BECs) with long-range anisotropic dipolar interactions 
and the interacting fermions in two-dimensions Q. 
Besides, the fundamental ideas of PTs have been utilized 
to study some models in quantum computation theory 
(QCT). Techniques based on PTs have also been appHed 
to engineer interesting Hamiltonians, whereby the out- 
standing perturbation theory gadgets (PTGs) technique 
was developed d, [7| ■ A short but incisive review of the 
history of PT and its new development in quantum many- 
body theory, quantum computation and quantum com- 
;^exity theory has been given recently by M. M. Wolf 

On the other hand, nonperturbation theories (NPTs) 
have been established to solve real scientific problems 
with the lack of a small parameter that acts as a pertur- 
bation parameter. In this case, some parameters will be 
artificially introduced as some formal perturbation pa- 
rameters and/or nonperturbation parameters. Usually, 
one of these artificial perturbation parameters have to 
be set finite, say, the unit one, while others will be deter- 
mined in some nonperturbative ways such as some types 
of optimized approaches. 

Considering the fact that physical quantities should be 



independent of any particular perturbation and/or non- 
perturbation method used to calculate them, Stevenson 
9] proposed an optimized perturbation theory (OPT) 
with the principle of minimal sensitivity (PMS) as a key 
ingredient. Some other perturbation methods appeared 
in 1970's and 1980's, such as the non-parameter expan- 
sion method [ ip| . and the linear 5 expansion method 
(LDE) [HI, [Illli, Hi, m, m, Ii3> have the same idea 
as implied in the OPT. These methods, especially the 
LDE, have found successful applications in many physi- 
cal contexts during the past three decades. 

Actually, as we will show later in this Letter, the ba- 
sic idea underlying those methods mentioned above can 
be integrated into a general framework of the (linear) 
homotopy analysis method (HAM) [l^ [l^ . In the lan- 
guage of field theory, it is possible to construct a new 
action Sq = T-C{So, S, q) by building a homotopy rela- 
tion between an original action S and a trial ideal action 
So (which is solvable and reflects as much as possible 
physics of S). Here S can represent any object or quan- 
tity for any scientific problem, such as the Green func- 
tion, density operator, distribution function, generating 
functional, differential equations, and so on. q G [0, 1] 
is called a homotopy parameter. Sq is fixed at the two 
endpoints, i.e., 5*^=0 = Sq and Sq^i = S, while can pos- 
sess different expressions when q ^ 0,1, depending on 
7^(5*0, S", g). Hence, different homotopy relations yield 
different actions. For example, a linear homotopy rela- 
tion yields an action as Sq = {1 — q)So + qS, which is 
exactly the form used in the LDE and other methods 
mentioned above. However, the homotopy relation is by 
no means constrained to be linear. We can introduce a 
simple and direct nonlinear extension 

n / ^ \ 

Sq = (l-<z)5o+^q^5,+<z"+M5-^5, .(1) 
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In the standard procedure of the LDE, 5*0 contains 
one or more auxiliary parameters rji. In our general- 
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ized nonlinear homotopy relation ([T]), auxiliary actions 
S'l, 5*2, . . . , S'n also contain some parameters These 
parameters are also artificial, which means that Sq=i — S 
is independent of both rji and . This can be easily ver- 
ified with the help of Eq. ([1]), as all 5*^ {i = 1 . . .n) will 
be eliminated when the nonlinear homotopy parameter 
q is set to 1. In order to describe it more explicit and 
without losing the nonlinear property of Eq. JT]), let us 
throw away 0{q^) terms to obtain 

Sg = {l-q)So{v)+qSi{0+q'{S-S,{0), (2) 

with only two auxiliary parameters rj and ^. Any desired 
physical quantity $ can be evaluated as a perturbation 
series of q, which is set equal to 1 at the final step. If 
the series is truncated at q^,N > n + l = 2, the ap- 
proximant is obviously dependent of rj and ^, while 
it will have no relation with those two auxiliary parame- 
ters when N — » oo. In practice, N is always finite, it is 
therefore necessary to choose these parameters nonper- 
turbatively to make the expansion reasonable. We will 
take the PMS as an effective criterion. According to the 
PMS, the auxiliary parameters 77 and f should be selected 
to minimize the sensitivity of the perturbation expansion 
to some small variations in them. It is remarkable that a 
nonlinear homotopy relation and the combination of the 
perturbation and nonperturbation theories via the PMS 
arc two fundamental ingredients of our nonscnsitivc non- 
linear homotopy approach (NNHA). 

Before moving to concrete examples, let us discuss a 
little more about the PMS. There is a natural way to fix 
the parameters ^ and 77, which is to solve the following 
equation system 



dr] 



= 0, 



c)$(^)(??,C) 
9? 



= 



(3) 



similar to the usual OPT theory. However, the above 
system does not always have one solution. Eqs. ([3]) might 
have many solutions {^i, ry^}, therefore, we have to find a 
method to fix one of them according to the PMS. To this 
end, we define the Af-th sensitivity k 
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(4) 



where 



denotes the norm. It is clear that the sensitiv- 



ities K^^'' for all M should tend to zero as ^ cxo, be- 
cause the real model is {^, ?/} independent. Consequently, 
it is natural to pick out the final required , 77} by min- 
imizing the first few sensitivities — 1,2,..., k) for 
the smallest k. More details can be seen in the following 
concrete examples. It is clear that a solution (if it exists) 
of Eqs. ([3]) is really a minimum of the first sensitivity 
(FS) ^[""l 

Now we use our theory to study an eigenvalue problem 
of a given Hamiltonian. This is a basic problem included 
in almost all quantum physical problems. For instance. 



various models in condensed matter physics deal with 
ground-states of spin Hamiltonians, and quantum com- 
putation can be performed by encoding a solution of a 
computation problem into the ground-state of a Hamil- 
tonian Q . The calculation of energy levels for a Hamilto- 
nian system is also a crucial problem in solid physics 0] ■ 
As the first example, we apply the theory to calculate 
the ground-state energy of a sextic quantum anharmonic 
oscillator 



(5) 



Generally, through a nonlinear homotopy relation in the 
form of Eq. we can construct a new Hamiltonian 



Hg^il- q){Ho) + q{Hi) +q\H- Hi). 



(6) 



Here we choose Hq ~ + lux^ and Hi = + + \x^ , 
and thus obtain 



H{q,UJ,\) 



(1 - q){p^ + UJX^) + q{p^ +x^ + \x^) 
+g2(l-A)x^ (7) 



It is obvious that Hq is the Hamiltonian of a harmonic 
oscillator with an auxiliary parameter uj denoting the fre- 
quency of the oscillation. Hi is a new term which can 
not be included without the second order term of q. It 
is an auxiliary Hamiltonian with a parameter A denoting 
the intensity of an additional sextic effect. It is easy to 
verify that H{l,Ljj,a) ~ H, and H{0,uj,a) — Hq. 

Making advantage of Bender and Wu's work we 
obtain the ground-state energy of H{q,uj,a) as a power 
series of q. 



(15A - 46.4 + 4.^2) 

EHT{q,(^,\) = ^ ^^^^3 q 



(8) 



which can be easily calculated to any order of q, with 
the help of the symbolic calculation softwares. Using 
MAPLE, we calculate the series up to order 21, and ob- 
tain an analytical expression of the approximate ground- 
state energy E^j,' {q,uj. A), and finally set q = 1. 

Now comes the problem to determine parameters lo 
and A under the PMS. Roughly speaking, the PMS 

tends to pick out parameters around which the profile 
f 21) 

of E)^:^'{1, UJ, A) is the "flattest" . So what can be learned 
from the profile of Ef^jl as displayed Fig. 1, where the 

(21) 

symbolic E denotes E^j^j!? It is discovered that when 
{uj,X) varies in the region [3.7,4] x [0.25,0.4], the varia- 
tion in the energy E^j^^ is less than 10~^. We take a very 
accurate value Eex = 1.435 624 619 as an "exact" result, 
which is given by Meifiner and Steinborn [2l[. Thus we 



have an error E 



(21) 
HT 



Eex vary from —4.62 x 10 



to 



5.38 X 10 ^ in the above region. In the case that the 
exact value of a solution is not known, we can calculate 
E^^ — E^j^^ as an error. It is reasonably found that 



E 



(21) 
HT 



E 



(20) 
HT 



10' 



To pick out an optimal point in 



this region, let us go further from the figure to determine 



FIG. 1: -E^y as a function of uj and A in the region [3.7, 4] x 
[0.25,0.4]. 
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FIG. 2: The curves of ^ q (sohd) and ^^fS- = 

(dash). 



a pair of nonsensitive parameters numerically. In Fig. 2 

we draw the curves of = (solid) and — gf^ = 

(dash), and find that they do not intersect in this region. 
In this situation, we can not obtain a pair of {uj, A) to 
vanish these two first-order partial derivatives simulta- 
neously. Motivated by the spirit of the PMS, we mini- 
mize the FS for = E%}^{uj, A) by simply taking the 
norm as = . Consequently, a pair of parameters 
u) = 3.847 166 and A = 0.377 585 is obtained to mini- 
mize the FS. The substitution of these parameters into 
the approximant (1, w. A) leads to an approximate 
energy E^^^ — 1.435 624 with an error of —5.53 x 10^''. 

As a compare, we also give an approximant truncated 
at order 21 for the ground-state energy by using a linear 
homotopy relation as iJ(q, w) — [l — q){p'^+LLix^) + q{p'^ + 

+ x^), which has been applied in the LDE method. It 

yields an approximant Ej^jJ^{q,uj) which only depends 
on w, after q is set to be 1. After some similar proce- 

(21) 

dures, we obtain a final result i?)^^^^(l, cjq) = 1-435 428 
with luq = 5.254 070 satisfying dE^'^^ /doj = 0. The error 
of E^f^^ is —1.97 X 10~^, which is about 10^ times larger 

(21) 

than that of E\jj! . The above example and the cor- 
responding numerical results show that our theory does 
largely improve the accuracy of the approximant without 
more difficulty. We have also study the ground and ex- 
cited states for some quartic and octic anharmonic oscil- 
lators in the same way, and drawn the same conclusions. 
It is remarkable that when the coefficients of the anhar- 
monic terms are large, the method is still applicable and 
the accuracy is at the same level. 

Next we go further to study non-hcrmitian Hamil- 
tonian systems. Recently, many progresses have been 



FIG. 3: S^y' as a function of ui and A in the region 
[1.68,1.69] X [2.153,2.175]. 

made for non-hermitian PT-symmetric physics. We are 
thus motivated to investigate the energy spectrum of a 
non-hermitian Hamiltonian with PT-symmetry. Bender 
and his co-workers have demonstrated that such a VT- 
symmetric Hamiltonian 

i?' = + + (9) 

has a real and positive spectrum [22]. According to our 
theory, the nonlinear homotopy relation can also be built 
to link a non-hermitian Hamiltonian with some auxil- 
iary Hamiltonians that can be either hermitian or non- 
hermitian. For instance, a nonlinear homotopy Hamilto- 
nian can be constructed as 

H'{q,LO,\) = [l-q)[p^ +LOx'^)+q{p^ + ^x'^ +i\x'^) 

+iq^{l-\)x^. (10) 

It can be seen that H'{q, w, A) also has the PT-symmetry, 
and its ground-state eigenvalue can be expressed as 

E'HT{q.^.^)^^+ g + (11) 

where the coefficients of each order of q are real. Truncat- 
ing the series at order 21 and setting q = 1 yield an ap- 
proximant E'I^^^ for the ground-state energy of H' . The 
pair that eliminates the FS (g]) for = E'^^\uj, A) is 

{uj ^ 2.160 887, A = 1.684 022}, and the profile of E^^^ 
around this stationary point is shown in Fig. 3 with E 
denoting E'j^^\ At the stationary point, E'^r^^ equals 
1.167 384, with an error 7.61 x 10~^ compared with the 
numerical result 1.167 46 given by Bender and Dunne 

This example reveals great validity of our theory on 
investigating non-hermitian "PT-symmetric Hamiltonian 
systems. Besides, our theory can be used to study more 
complicated problems beyond a linear eigenvalue prob- 
lem. It is shown in a series of recent works that non- 
linear localized structures can be generated in VT pe- 
riodic potentials [13, |2^ . Optical solitons are produced 
by adding a self-focusing Kerr nonlinear PT-symmetric 
potential iW{x), which satisfies W{—x) = —W{x), to 
the usual nonlinear Schrodinger equation. Our theory 
can also be applied to investigate this type of nonlinear 



4 



PT-symmetric system by building suitable homotopy re- 
lations associating a system hard to be exactly solved 
with those solvable. In this way, the knowledge of the 
latter can be used to explore the original system approx- 
imately, and the necessary accuracy is entirely ensured 
by the principle of minimal sensitivity. 

As is known, physical models that can be solved ex- 
actly are very rare. In a recent work on the structure 
of trapped Bose- Einstein condensates (BECs) with long- 
range anisotropic dipolar interactions, it is found that a 
small perturbation in the trapping potential can lead to 
dramatic changes in the condensate's density profile for 
sufficiently large dipolar interaction strengths and trap 
aspect ratios The method the authors used is the tra- 
ditional perturbation method which require the pertur- 
bation to be small, our theory can thus be used to extend 
their work to the range where the perturbation is larger. 
So matter-wave systems with dipolar interactions would 
be a suitable field for our theory. Besides, the theory can 
also facilitate the study of quantum many-body systems, 
where perturbation techniques are essential tools [8|. 

To conclude, we have proposed a novel approach called 
NNHA, a combination theory of perturbation and non- 
perturbation with the help of the PMS and nonlinear 
homotopy realization. The new method relies on a non- 



linear homotopy to construct an approximant with aux- 
iliary parameters that does not exist in the original mod- 
els. The nonlinear homotopy parameter q is taken as 
a perturbation parameter and finally fixed to 1 as usual, 
because the original model is related to q = 1 only. Other 
parameters are nonperturbatively fixed at the end of the 
calculation by implementing the PMS, which requires the 
approximant have the least dependence on these auxil- 
iary parameters for the reason that the original system 
is independent of them. The theory has been applied 
to study the energy spectrum of some hermitian Hamil- 
tonian systems as well as non-hermitian PT-symmetric 
Hamiltonians. Highly accurate numerical results demon- 
strate the validity of the theory. Possible applications 
can be made further to PT-symmetric optical systems, 
dipolar Bose- Einstein condensates (BECs) , and quantum 
many-body systems. Actually, we do believe that the 
method can be used to solve any scientific problem where 
mathematics is required. 
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